Renormalization group approach to 2D Coulomb interacting Dirac fermions with 

random gauge potential 
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We argue that massless Dirac particles in two spatial dimensions with 1/r Coulomb repulsion 
and quenched random gauge field are described by a manifold of fixed points which can be accessed 
perturbatively in disorder and interaction strength, thereby confirming and extending the results 
of arXiv:0707.4171 At small interaction and small randomness, there is an infra-red stable fixed 
curve which merges with the strongly interacting infra-red unstable line at a critical endpoint, along 
which the dynamical critical exponent z = 1. 



The properties of two dimensional massless Dirac 
fermions have recently sprung back into focus, largely due 
to the experimental discovery of the quantum Hall effect 
in graphene[l|,0], the single layer graphite. Moreover, the 
ability to control the density of carriers by the electrical 
field effect allows experimental access to the rich physics 
of the neutrality point, where in the clean non- interacting 
picture, the conduction and the valence bands touch. It is 
well known Q that, at the neutrality point, the exchange 
self-energy gives a logarithmic enhancement of the Fermi 
velocity vp ^ vp + (e^/4ed) ln(A/fc) where fc is a small 
wavevector near the nodal point's] and is the dielectric 
constant of the medium. Physically, this effect is due to 
the lack of screening of the 1/r Coulomb interaction, an 
important consequence of which is the suppression of the 
single particle density of states {N{E)) at low energies. 
This in turn leads to the suppression of the electronic 
contribution to the low temperature specific heatQ. 

This suppression of N{E) may lure one into the (in- 
correct) conclusion that, at T = 0, the Coulomb inter- 
actions turn the clean system into an electrical insula- 
tor. However, the vertex corrections contribute an ex- 
actly compensating enhancement of the conductivity 
making the system a metal with its residual conduc- 
tivity asymptotically equal to the non-interacting value 
CTo = {n/8)e^/h per node. 

In this work, wc analyze the effects of the unscreened 
Coulomb interactions and the quenched random gauge 
disorder beyond leading order in the perturbative renor- 
malization group (RG)ofRef. d. Our principle findings, 
which support and extend those of Ref. Q are twofold: 
first, in the clean case, there is an unstable fixed point 
at finite strength of Coulomb interactions characterized 
by the dimensionless ratio a = e'^ / (hedVp) which rep- 
resents a quantum critical point (QCP) separating the 
semimetal from an excitonic insulator; and second, the 
interplay between Coulomb interactions and disorder in- 
duces a downward curvature of the fixed line [5, 6], caus- 
ing it to end at the clean QCP (see Fig.l). 

In two dimensions, the Hamiltonian for Coulomb in- 
teracting massless Dirac fermions in the presence of a 



FIG. 1: The renormalization group fiow diagram in the (dis- 
order) A - (interaction) a = e'^ /(edVp) plane. There is a line 
of stable fixed points at small A and small a which merges 
with the line of unstable fixed points at the critical end point. 
The (clean) unstable fixed point at Oc corresponds to a quan- 
tum phase transition into an excitonic insulator. Above the 
critical A*, the disordered but non-interacting fixed line is 
unstable directly to the insulator. 



quenched random gauge field can be written as 

n = Ho + + V (1) 

with the free part given by, 

N 



^o = E / ^'"^"^ [V'l(r)p-aV'j(r) 



(2) 



where the operator ijjj annihilates a two-component Dirac 
fermion, p = — iV, Ui is the ith Pauli matrix, and we 
have set h = ks = 1 for convenience, unless otherwise 
stated. TV represents the number of fermion species; for 
single layer graphene, = 4. The disorder part of the 
Hamiltonian is 



N 



J2 ^j^F I d'r 



a^(r)^ (r)a^V,(r) (3) 



where /i = 1 , 2 andthe quenched random gauge field is 
assumed to be uncorrelated 



The charges Xj = (—1)-' vary in sign from node to node as 
dictated by the overall time reversal symmetry of the sys- 
tem. The Coulomb interaction between Dirac fermions 
is given by 



V = 



1 



d^rcfr' 



5n{r) 



r — r' 



-dh{r') 



(5) 



where 6n{r) = J2j V'j(r)'(/'i(r) — "-o and e is the electronic 
charge. The background charge density no ensures over- 
all charge neutrality. 

This system can be described by three "coupling con- 
stants" - e^, A and vp - oi which A is dimensionless 
and, in our units, both vp and e^/e^ have dimensions of 
velocity. In what follows, we will set ed — 1 and restore it 
in the final results by rescaling the charge. Generically, 
these coupling constants flow under the renormalization 
group transformation. However, the charge does not 
flow because it is a coefficient on a non-analytic term 
in the action, and as detailed elsewhere A doesn't flow 
either The entire flow of the renormalized coupling 
constants then comes from the scale dependence of the 
Fermi velocity. The RG beta functions take the form: 
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(6) 
(7) 
.(8) 



where the ellipses mean terms of cubic order in the double 
expansion in small and A. The lowest order terms 
in the expansion come from the exchange diagrams for 
disorder and interactions (see Fig. 2), respectively [HI . 
The result of our analysis presented below is the following 
values of the above coefficients 
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The vanishing coefficient A agrees with the result of Lud- 
wig ei.a/.Q that, for = 0, the dynamical critical expo- 
nent z = 1 — dlnvp/dliiK — 1 — A/tt holds to all orders 
in perturbation theory. The values of B and C are the 
new results first reported in this paper. 

If we rescale hy vp, we can define a dimensionless 
coupling constant ap — /[hedVp) which characterizes 
the strength of Coulomb interactions. The corresponding 
flow diagram is shown in Fig.l. In the clean case, small 
ap flows to zero due to the growth of the Fermi velocity. 
Ki ap = ac, there is a quantum phase transition into an 
excitonic insulator, controlled by the strongly interacting 



FIG. 2: Diagrams contributing to the self energy to the 
f'order in ap and Aa- As can be readily seen from the link 
cluster expansion, the second diagram comes with an overall 
- sign. 




FIG. 3: Diagrams contributing to the self energy to the order 
a% (first row) and (second row). 



flxed point. Within the above approximation, for iV = 4, 
which is appropriate for the single layer graphenc, the 
unstable fixed point appears at ac = 1/4C « 0.833. As 
discussed in greater detail below, since the expansion of 
the beta function ([8]) is carried out in apN ^ 1, the 
semimetal-insulator fixed unstable point appears beyond 
the strict validity of the perturbative RG. Nevertheless, 
for finite A, there is a fixed manifold given by 



^Cap) 



4 l + TiBap 



(10) 



which is shown in Fig.l by the solid (red) curve. Impor- 
tantly, at small ap and A, this manifold represents the 
line of stable fixed points p, HI which is asymptotically 
exacts and which merges with the line of unstable fixed 
points at the critical endpoint (a*. A*). 

Next, we detail the perturbative renormalization group 
calculation which we perform in dimensional regulariza- 
tion scheme by analytically continuing the space integrals 
io D — 2 — e where e > 0. This formal device serves 
as a regulator for various divergent integrals. The bare 
Green's function has the form lal 



Go{iuj, k) = {—iiv -I- cr • k) ^ ^ 



iuj + a - li 
^2 + k2 



The resulting self-energy matrix, which is defined 
through the two point irreducible vertex T = Q^^ as 



Tk{iLu) = —ibj + vpa ■ k + I]k(«^), 



(11) 



is calculated at finite external momentum k and fre- 
quency iuj. To illustrate the procedure, consider the self- 
energy to first order in the coupling constant (the first 





FIG. 4: Diagrams contributing to the self energy to the order 



diagram in Figl2]): 



r[i-f]r[i2±i]r[i^] 



{An)' 



r[D] 



■ kfc-^ + ^ ln(647re^^^)^ . (12) 

In the last line, we have used the analytic properties 
of the Euler F— functions and expanded the result near 
e = 0. The Euler constant je « 0.577. The pole at 
D = 2 corresponds to a logarithmic divergence when a 
large momentum cutoff is introduced directly in two di- 
mensions i.e. the Hartree-Fock approximation leads to 
the logarithmic enhancement of velocity and Sk(iti-') — 
■ kin A/A: is frequency independent. This logarith- 
mic enhancement of the Dirac velocity persists to second 
order in the coupliiig constant expansion, in qualitative 
agreement with Ref|9|, as well as the leading order expan- 
sion in large N [loj . 



Similarly, the leading order self energy due to the dis- 
order scattering (second diagram in Fig. [2]) is 
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13) 



Again, the pole aX D = 2 corresponds to a logarithmic 
divergence of the self-energy. 

Before addressing the higher order contributions to the 
self-energy, let us define the renormalization conditions. 
The standard relationship between the renormalized two 
point function rj^'^"(icj) and the bare one PT|) is 



(14) 



where Z is the wavefunction renormalization [11]. The 
renormalized coupling constants can now be defined 
through the following renormalization conditions 
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(15) 



(16) 



Physically the above equations demand that at the 
renormalization scale k, the renormalized single particle 
Green's function G'^'^{iu)) takes the form 



^ren - 



From the above considerations, it is clear that, to find the RG flows, we need only to find the self energies for 
= |k| = K. At the orders A^, e^A, and we find 
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(17) 
(18) 
(19) 



Thus, to this order in the perturbative expansion, the 
self energy matrix in Eq. (fTT|) is 

Ek(«w) = S^^'='*(^w)^-I]^^"(^t^)-HI]£(^^^)-KS^(^t^)-t-S^i(^CJ), 



and using Eqs. (11411161) we find the RG flow equation for 
the renormalized Fermi velocity to be 
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This is the result displayed in Eq 
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It is apparent from 



this flow equation that at A = the Fermi velocity in- 
creases logarithmically provided that the dimensionless 
coupling a < Uc- Such logarithmic growth implies sup- 
pression of the electronic density of states near the Dirac 
point and concomitant suppression of the speciflc heatQ. 

In the case of a > ac, the runaway RG flows can be in- 
terpreted as the flow towards an excitonic insulator. The 
physical nature of this insulator depends on the details 
of the lattice model. For spinless fermions on a hon- 
eycomb lattice, for example, it would correspond to a 



state with a spontaneous breaking of (plaquette or bond 
centered) inversion symmetry and unequal population of 
the two sublattices; for strictly 1/r Coulomb repulsion 
and for the spinfuU case it would correspond to an " an- 
tiferromagnetic" order (but no unit cell doubling) with 
unequal spin population of the two sublattices, although 
details of the short range part of the repulsive interac- 
tions could destabilize it towards the inversion symmetry 
broken state [l^. 

We conclude with a discussion of the validity of the 
approximations employed here. First, the RG was or- 
ganized perturbatively in both a and A so the weak 
coupling portion of the fixed line in Fig[l] is rigorously 
justified. Additionally, the downward curvature of this 
line implied by the sign of the terms in the second or- 
der expansion is also rigorous. Diagramatically, inte- 
gration of these RG equations corresponds to an infinite 
parquet-like resummation of the two leading logarithms 
at each order in the pertrubative expansion (see for in- 
stance [l3|). 

On the other hand, the unstable portion of the fixed 
line and the (unstable) clean fixed point appear at a 
finite value of the coupling constant at which apN is 
not small. They are, therefore, beyond the reach of the 
perturbative RG. Nevertheless, the very existence of the 
clean unstable fixed point is perhaps on somewhat firmer 
footing 14, 15, 1^, and due to the negative curvature of 



the IR stable fixed line in the perturbatively accessible 
region of the fiow diagram, we expect that the gross topo- 
logical features of the fixed manifold in FiglU if not its 
quantitative aspects, are valid. 

We wish to thank Professors I. Herbut and Z. 
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reading of the manuscript. 



[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 

A. A. Firsov, Nature 438, 197 (2005). 
[2] Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 

438, 201 (2005). 
[3] J. Gonzales, F. Guinea, and M. A. H. Vozmediano, Nucl. 

Phys. B 424, 595 (1994). 
[4] O. Vafek, Physical Review Letters 98, 216401 (pages 4) 

(2007). 

[5] I. F. Herbut, V. .Juricic, and O. Vafek (2007), URL 
http : / /www . arxiv . org/abs/0707 .4171 

[6] T. Stauber, F. Guinea, and M. A. H. Vozmediano, Phys- 
ical Review B 71, 041406 (2005). 

[7] A. Ludwig, M. P. A. Fisher, R. Shankar, and G. Grin- 
stein, Phys. Rev. B 50, 7526 (1994). 

[8] J. W. Negele and H. Orland, Quantum Many-Parttcle 
Systems (Addison- Wesley, Reading, MA, USA, 1988). 

[9] E. G. Mishchenko, Physical Review Letters 98, 216801 
(2007). 

[10] J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Phys. 

Rev. B 59, R2474 (1999). 
[11] M. E. Peskin and D. V. Schroeder, An Introduction to 

quantum field theory (Addison- Wesley, Reading, MA, 

USA, 1995). 

[12] I. F. Herbut, Physical Review Letters 97, 146401 (2006). 
[13] I. F. Herbut, A Modem Approach To Critical Phenomena 

(Cambridge University Press, Cambridge, UK, 2007). 
[14] D. V. Khveshchenko, Phys. Rev. Lett. 87, 246802 (2001). 
[15] E. V. Gorbar, V. P. Gusynin, V. A. Miransky, and I. A. 

Shovkovy, Phys. Rev. B 66, 045108 (2002). 
[16] H. Leal and D. V. Khveshchenko, Nucl. Phys. B 687, 323 

(2004). 



